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Abstract 

A library for reading and writing data in the SUSY Les Houches Accord 2 format 

■ is presented. The implementation is in native Fortran 77. The data are contained in 
. a single array conveniently indexed by preprocessor statements. 

o ■ 
in : 

\o '• 1 Introduction 

o ; 

The original SUSY Les Houches Accord [1 (SLHA1 in the following) has standardized and 
significantly simplified the exchange of MSSM input and output parameters between such 
disparate applications as spectrum calculators and event generators. Meanwhile, agreement 

■ has been reached also about the encoding of many extensions of the MSSM which has led 
>• . to a preliminary SLHA2 document [2]. 

^ . While the SLHA specifications include the precise formats for Fortran I/O, it is nev- 

ertheless not entirely straightforward to read or write a file in SLHA format. The present 
library provides the user with simple routines to read and write files in SLHA format, as 
well as a few utility routines. One thing the library does not do is modify the numbers, 
which means there is no routine to compute, say, a particular quantity at a new scale. The 
data structures and subroutines are set up such that only very few changes are necessary 
when upgrading from the SLHALib 1 [3] . 

Sect. 121 describes the organization of the data structures, Sect. E3 gives the reference 
information for the library routines, Sect. 0] shows the usage in some examples, Sect. El 
contains download and build instructions, and Sect. |H1 summarizes. 



2 Data structures 

The SLHA library is written in Fortran 77. All routines operate on a double complex 
array, slhadata, which is about the simplest conceivable data format for this purpose in 
Fortran. For convenience of use, this array is accessed via preprocessor statements, so the 
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user never needs to memorize any actual indices for the slhadata array. A file containing 
the preprocessor definitions must thus be included. 

The slhadata array consists of a 'static' part containing the information from SLHA 
BLOCK sections and a 'dynamic' part containing the information from SLHA DECAY sections. 
The static part is indexed by preprocessor variables defined in SLHA . h, the dynamic part is 
accessed through the SLHANewDecay, SLHAFindDecay, SLHAAddDecay, SLHAGetDecay, and 
SLHADecayTable functions and subroutines (see Sect. |3J). 

In addition, descriptive names for the PDG codes of the particles are declared in PDG . h. 
These are needed e.g. to access the decay information. 

2.1 SLHA blocks 

The explicit indexing of the slhadata need not (and should not) be done by the user. 
Rather, the members of the SLHA data structure are accessed through preprocessor vari- 
ables. Tables ITtil 01 list the preprocessor variables defined in SLHA . h which follow closely the 
definition of the Accord. Note that preprocessor symbols are case sensitive. On the down- 
side, there is no way to guard against out-of-range indices, not even with compiler flags. 
This is because the preprocessor has no such checks and the compiler cannot determine a 
posteriori whether the single index it sees addresses the 'right' part of the array. 

As far as there is overlap, the names for the block members have been chosen similar 
to the ones used in the MSSM model file of FeynArts 0]. Following is a list of common 
index conventions. This is only for a rough orientation: the actual indices and their ranges 
are always given explicitly in the Tables. 

t — 1...4 (s)fermion type: 1 = (s)neutrinos, 

2 = isospin-down (s)leptons, 

3 = isospin-up (s) quarks, 

4 = isospin-down (s) quarks 

(s)fermion generation 

number of sfermion mass-eigenstate, 

in the absence of mixing 1 = L, 2 = R 

c = 1 . . . 2 number of chargino mass-eigenstate 

n = 1 ... 4 number of neutralino mass-eigenstate 

For each block B the offset into slhadata and the length are respectively defined as 
Off set5 and LengthS. The contents of the block can be addressed through the macro 
BlockS (z), where i runs from 1 to LengthS. 

Matrices have a "Flat" array superimposed for convenience, in Fortran's stan- 
dard column-major convention, e.g. USf (1,1) = USfFlat(l), USf (2,1) = USfFlat(2), 
USf (1,2) = USf Flat (3), USf (2,2) = USf Flat (4). This makes it possible to e.g. copy 
such a matrix with just a single do-loop. 
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Block name 


Array and length 


Members 


oJ^lxNr (J 


DlOCKorllllO \Ti) 


SPInf o_Severity 






Lengxnor ini o 


SPInf o_NLines 








SPInf o_Code(n) 


n = 1 ... 15 






oJrllll 0_ I GXX \ l } 7iJ 


qo nor aaoress airecriy 


DCTNFO 


BlockDCInf o (n) 


DCTnf o Spvpti tv 

LJ Kj _L XJ L W h_;GVCL _L L> y 






LengthDCInf o 


DCTnf o NT inpq 








DCTnfo CorfpCr)") 

J-/ \J X XX J. W \j\J\JLKZ \ I v J 


n = 1 ... 15 






DPTrifn Tpirt (i n 1 

i/vi. ±1 J. VJ 1 C A U \i j / ^ / 


Ho not prlrlrptit; Hirprtlv 


MUL'orjL 


b±oc.ki v iucloG-L 


ModSel_Model 






J_i G IlgX LLrl Qo G X 


ModSel_Content 








ModSel_RPV 








hit in n /ttit t 

ModSel_CPV 








Tiff _ JO— .~\ T - 'T T 

ModSel_FV 








Tiff TO TO* ITij. 

ModSel_GridPts 








ModSel_Qmax 








ModSel_PDG(i) 


i = 1 ... 5 


SMINPUTS 


BlockSMInputs(n) 


SMInputs_AlfaMZ 






LengthSMInputs 


SMInputs_invAlf aMZ 


ee SMInputs_AlfaMZ 






SMInputs_GF 








SMInputs_AlfasMZ 








SMInputs_MZ 








SMInputs_Mf (t,g) 


t = 2...4, g=1...3 






SMInputs_Mf Flat ( i) 


i = 1...9 






SMInputs_Me 


= SMInputs_Mf (2,1) 






SMInputs_Mu 


= SMInputs_Mf (3,1) 






SMInputs_Md 


= SMInputs_Mf (4,1) 






SMInputs_Mmu 


= SMInputs_Mf (2,2) 






SMInputs_Mc 


= SMInputs_Mf (3,2) 






SMInputs_Ms 


= SMInputs_Mf (4,2) 






SMInputs_Mtau 


= SMInputs_Mf (2,3) 






SMInputs_Mt 


= SMInputs_Mf (3,3) 






SMInputs_Mb 


= SMInputs_Mf (4,3) 



Table 1: Preprocessor variables defined in SLHA.h to access the slhadata array. The 
equivalence symbol (=) indicates that the l.h.s. is just an alias for the r.h.s., not a new 
variable. 
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Block name 


Array and length 


Members 


MINPAR 


BlockMinPar (n) 


MinPar_Q 








L e ngt hM i nP ar 


MinPar_M0 










MinPar_Lambda 


= 


MinPar_M0 






MinPar_M12 










MinPar_Mmess 


= 


MinPar_M12 






MinPar_M32 


= 


MinPar_M12 






MinPar_TB 










MinPar_signMUE 










MinPar_A 










MinPar_N5 




MinPar_A 






MinPar_cgrav 






EXTPAR 


BlockExtPar(n) 


ExtPar_Q 








LengthExtPar 


ExtPar_Ml 










ExtPar_M2 










ExtPar_M3 










ExtPar_Af (t) 


t = 


= 2. ..4 






ExtPar_Atau 


= 


ExtPar_Af (2) 






ExtPar_At 


= 


ExtPar_Af (3) 






ExtPar_Ab 


= 


ExtPar_Af (4) 






ExtPar_MHu2 










ExtPar_MHd2 










ExtPar_MUE 










ExtPar_MA02 










ExtPar_TB 










ExtPar_MSS(#,<7) 


g - 


= 1...3, q= 1...5 






ExtPar_MSL(g) 




ExtPar_MSS(5f,l) 






ExtPar_MSE(#) 




ExtPar_MSS(g,2) 






ExtPar_MSQ(gO 




ExtPar_MSS(g,3) 






ExtPar_MSU(#) 




ExtPar_MSS(5,4) 






ExtPar_MSD(#) 




ExtPar_MSS(5,5) 






ExtPar_N5(#) 


g -- 


= 1...3 






ExtPar_lambda 










ExtPar_kappa 










ExtPar_Alambda 










ExtPar_Akappa 










ExtPar_MUEeff 







Table 2: Preprocessor variables defined in SLHA.h to access the slhadata array (cont'd). 



4 



Block name 


Array and length 


Members 


MASS 


BlockMass(n) 


Mass_Mf (t ,g) 


t=l... 


4, o = l. 


.3 




LengthMass 


Mass_MfFlat(z) 


i = 1... 


12 








Mass_MSf (s,t,g) 


8 = 1.. 


2, t = 1. 


.4, 5 = 1. ..3 






Mass_MSfFlat(i) 


i = 1... 


24 








Mass_MZ 












Mass_MW 












Mass_MhO 












Mass_MHH 












Mass_MAO 












Mass_MHp 












Mass_MHl 


= Mass 


_MhO 








Mass_MH2 


= Mass 


_MHH 








Mass_MH3 












Mass_MAl 


= Mass 


_MA0 








Mass_MA2 












Mass_MNeu(n) 


n= 1.. 


.5 








Mass_MCha(c) 


c = l.. 


2 








Mass_MGl 












Mass_MGrav 








NMIX 


BlockNMix(n) 


NMix_ZNeu(ni ,n 2 ) 


ni,n 2 = 


1...4 






LengthNMix 


NMix_ZNeuFlat (i) 


i = l... 


16 




UMIX 


BlockUMix(n) 


UMix_UCha(ci,c 2 ) 


Cl , C 2 = 


1...2 






LengthUMix 


UMix_UChaFlat (i) 


i = 1... 


4 




VMIX 


BlockVMix(ra) 


VMix_VCha(ci,c 2 ) 


Cl,C 2 = 


1...2 






LengthVMix 


VMix_VChaFlat (i) 


i = 1... 


4 





Table 3: Preprocessor variables defined in SLHA.h to access the slhadata array (cont'd). 
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Block name 


Array and length 


Members 




R1 nrV^f Mi y (n ~) 

DIULlvOllJlA \ 1 1 J 


Sf Mix_USf (s\ ,S2,t) 


Sl 


S2 = 1 • • • 2, t = 2 . . .4 




T ono-t-Vi^-f Mi y 


Sf Mi y TTSfFlatf? f} 


i — 

v — 


= 1 ... 4, t = 2 . . . 4 


STAUMIX 


BlockStauMix(n) 


StanMiY TTSf f si <?o") 

kj U CLLL1 1XA Uul V^l ) ^ 2i ' 




SfMix USfCsi so 2") 




LengthStauMix 


StauMix_USf Flat (i) 





Sf Mix_USf Flat ( i , 2) 


STOPMIX 


BlockStopMix(n) 


StopMix_USf (s\ , S2) 





SfMix_USf (si,S2,3) 




LengthStopMix 


StopMix_USf Flat (i) 


— 


SfMix_USfFlat(i,3) 


SBOTMIX 


BlockSbotMix(n) 


SbotMix_USf (s 1 ,s 2 ) 


— 


SfMix_USf (si,s 2 ,4) 




LengthSbotMix 


SbotMix_USfFlat(i) 


— 


SfMix_USfFlat(i,4) 


ALPHA 


BlockAlpha(n) 


ii-L UIlcl_ii± pilcL 








LengthAlpha 








HMIX 


rslOCKrinlX KJl) 


HMix_Q 








Lengtnnnix 


HMix_MUE 










HMix_TB 










HMix_VEV 










HMix_MA02 






GAUGE 


BlocktGauge (n) 


Gauge_Q 








LengthGauge 


Gauge_gl 










Gauge_g2 










Gauge_g3 






MSOFT 


BlockMSoft(ra) 


MSoft_Q 








LengthMSoft 


MSoft_Ml 










MSoft_M2 










MSoft_M3 










MSoft_MHu2 










MSoft_MHd2 










MSoft_MSS(#,g) 


9 ' 


= 1...3, q= 1...5 






MSoft_MSL(gf) 




MSoft_MSS(5,l) 






MSoft_MSE(5) 




MSoft_MSS(5,2) 






MSoft_MSQ(c/) 




MSoft_MSS(5,3) 






MSoft_MSU(5) 




MSoft_MSS(£?,4) 






MSoft_MSD(c/) 




MSoft_MSS(^,5) 



Table 4: Preprocessor variables defined in SLHA.h to access the slhadata array (cont'd). 
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Block name 


Array and length 


Members 




BlockAf (n) 

T on rr+" Vi A "F 
J_.eilgT. Ilil I 


Af_Q(r) t = 2...4 

If Ifl'n, /r„ ^ n, n n — 1 ^ + — 9 A 
AI_AI (,<7i ,^2 >D <7lj<72 — l.-.O, E — Z...4 

Af_AfFlat(i,t) i = 1 ... 9, c = 2...4 


AE 


BlockAe(n) 

T p-n cr1~ Vi A p 


Ae_Q = Af_Q(2) 
Ap AfCfli no") = Af Af foi no 9") 
Ae_AfFlat(i) = Af _Af Flat (i , 2) 
Ae_Atau = Ae_Af(3,3) 


AU 


BlockAu(n) 
T .pti P"l~h An 


Au_Q = Af_Q(3) 
Au Affoi oo) = Af AfCfli r/o 3) 
Au_AfFlat(i) = Af_AfFlat(i,3) 
Au_At = Au_Af(3,3) 


AD 


BlockAd(n) 
LenrthAd 


Ad_Q = Af_Q(4) 

Ad Afffli flol = Af AfC(7i On 4} 

Ad_AfFlat(i) = Af_AfFlat(i,4) 
Ad_Ab = Ad_Af(3,3) 




BlockYf (n) 

T rr-l-Vi V-F 

LGngxn i x 


Yf_Q(r) t = 2..A 

V-F h-f C si rt +\ rt rt 1 Q -f A 

Yi_Ai{gi,g2,t) <7i,<72 — l...o, t — _...4 
Yf_AfFlat(i,r) i = 1...9, t = 2. ..4 


YE 


BlockYe(n) 
LengthYe 


Ye_Q = Yf_Q(2) 
Ye_Yf( 5 i )52 ) = Yf_Yf (gi,g 2 ,2) 
Ye_YfFlat(i) = Yf _Yf Flat (i , 2) 
Ye_Atau = Ye_Yf(3,3) 


YU 


BlockYu(n) 
LengthYu 


Yu_Q = Yf_Q(3) 
Yu_Yf( 5 i, 52 ) = Yf_Yf(y-i,y- 2 ,3) 
Yu_YfFlat(i) = Yf_YfFlat(i,3) 
Yu_At = Yu_Yf(3,3) 


YD 


BlockYd(n) 
LengthYd 


Yd_Q = Yf_Q(4) 
Yd_Yf( 5l)52 ) = Yf_Yf (51,52,4) 
Yd_YfFlat(i) = Yf_YfFlat(i,4) 
Yd_Ab = Yd_Yf(3,3) 



Table 5: Preprocessor variables defined in SLHA.h to access the slhadata array (cont'd). 
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T^lnck TiaiTip 

1 J lw V IV llCllllVj 


A rrav flnd lpnp"t,Vi 


TVTptv. ISprs 


RVLAMBDAIN 


BlockRVLambdaln(n) 
LengthRVLambdaln 


RVLambdaIn_lambda (i , j ,k) i,j,k=1...3 
RVLambdaIn_lambdaFlat (i) i = 1 ... 27 


DAVT A AID T\ A 

K V L A M Jd D A 


DXOCKrtvLianiDaa v nj 
LengthRVLambda 


n, v j_i am d q a _ 

RVLambda_lambda(i,j ,k) i, j, k = 1 . . . 3 
RVLambda lambdaFlat 0") i = 1 27 


RVLAMBDAPIN 


BlockRVLambdaPIn(n) 
LengthRVLambdaPIn 


RVLambdaPIn_lambdaP(i,j,fc) i,j,k= 1...3 
RVLambdaPIn_lambdaPFlat (i) i = 1 ... 27 


K V L A M Jd D A P 


DXOCKrtvLianiDaar Vii/ 
LengthRVLambdaP 


n, V L. cLIIl D Q a. r _ y 

RVLambdaP_lambdaP(z, j , k) z, j, fc = 1 . . . 3 
RVLambdaP lambdaPFlat (?) 7 = 1 27 


RVLAMBDAPPIN 


BlockRVLambdaPPIn(n) 
LengthRVLambdaPPIn 


RVLambdaPPIn_lambdaPP(i, j,/0 i,j,fe = 1...3 
RVLambdaPPIn_lambdaPPFlat (i) i = 1 ... 27 


TAA7T A A /TTZ» T~\ A I » I > 

KVLAMBDArr 


DlOCKKvLaHDQarr v^/ 1 

LengthRVLambdaPP 


KvLamDaarr _y 

RVLambdaPP_lambdaPP(i, i,j,k= 1...3 
RVLambdaPP lambdaPPFlat (?'") ?' = 1 27 


RVAIN 


BlockRVAIn(n) 
LengthRVAIn 


RVAIn_A(i,j,fc) i,j,k = l...S 
RVAIn_AFlat (i) z = 1 ... 27 


RVA 


bIockKVAI. ri; 
LengthRVA 


T>\7 A n 

KVA_LJ 

RVA_A(i,j,fc) i,j,fc = 1...3 
RVA_AFlat(i) i = 1 ... 27 


RVAPIN 


BlockRVAPIn(n) 
LengthRVAPIn 


RVAPIn_AP(i,j,A;) fc = 1 . . . 3 
RVAPIn_APFlat (i) i = 1 ... 27 


RVAP 


BlockRVAP(n) 
LengthRVAP 


RVAP_Q 

RVAP_AP(i,j,fc) i,j,k=l...3 
RVAP_APFlat («) i = 1 ... 27 


RVAPPIN 


BlockRVAPPIn(n) 
LengthRVAPPIn 


RVAPPIn_APP(«,j,fe) i,j,fc = 1...3 
RVAPPIn_APPFlat ( i) i = 1 ... 27 


RVAPP 


BlockRVAPP 
LengthRVAPP 


RVAPP_Q 

RVAPP_APP(i,j,fc) z, j, = 1 ... 3 
RVAPP_APPFlat («) i = 1 ... 27 



Table 6: Preprocessor variables defined in SLHA.h to access the slhadata array (cont'd). 
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Block name 


Array and length 


Members 


RVKAPPAIN 


BlockRVKappaln(n) 
LengthRVKappaln 


RVKappaIn_kappa(i) i = 1 ... 3 


RVKAPPA 


BlockRVKappa(n) 
LengthRVKappa 


RVKappa_Q 

RVKappa_kappa(i) i = 1 ... 3 


RVDIN 


BlockRVDIn(n) 
LengthRVDIn 


RVDIn_D(i) i=1...3 


RVD 


BlockRVD(n) 
LengthRVD 


RVD_Q 

RVD_D(i) i=l... 3 


RVSNVEVIN 


BlockRVSnVEVIn(n) 
LengthRVSnVEVIn 


RVSnVEVIn_VEV (i) i = 1 . . . 3 


RVSNVEV 


BlockRVSnVEV(n) 
LengthRVSnVEV 


RVSnVEV_Q 

RVSnVEV_VEV(i) i=1...3 


RVMLHlSQIN 


BlockRVMLHlSqln(n) 
LengthRVMLHISqIn 


RVMLHlSqIn_MLH12 (i) i = 1 . . . 3 


RVMLHlSQ 


BlockRVMLHlSq(n) 
LengthRVMLHISq 


RVMLHlSq_Q 

RVMLHlSq_MLH12 (i) i = 1 . . . 3 


RVNMIX 


BlockRVNMix(n) 
LengthRVNMix 


RVNMix_ZNeu(ni,n 2 ) ni,n 2 = 1...7 
RVNMix_ZNeuFlat (i) i = 1 ... 49 


RVUMIX 


BlockRVUMix(n) 
LengthRVUMix 


RVUMix_UCha(ci,c 2 ) ci,c 2 = 1...5 
RVUMix_UChaFlat (i) i = 1 ... 25 


RVVMIX 


BlockRVVMix(n) 
LengthRVVMix 


RVVMix_VCha(ci,c 2 ) ci,c 2 = 1...5 
RVVMix_VChaFlat (i) i = 1 ... 25 


RVHMIX 


BlockRVHMix(n) 
LengthRVHMix 


RVUMix_UH(ft,i,/i 2 ) hi,h 2 = 1...5 
RVUMix_UHFlat (i) i = 1 ... 25 


RVAMIX 


BlockRVAMix (n) 
LengthRVAMix 


RVAMix_UA(/ii,/i 2 ) /ii,/i 2 = 1...5 
RVAMix_UAFlat (i) i = 1 ... 25 


RVLMIX 


BlockRVLMix (n) 
LengthRVLMix 


RVLMix_CLep(Zi,Z 2 ) /i,/ 2 = 1...8 
RVLMix_CLepFlat (i) i = 1 ... 64 



Table 7: Preprocessor variables defined in SLHA.h to access the slhadata array (cont 
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Block name 


Array and length 


Members 


VCKMINPUTS 


BlockVCKMInputs (n) 
LengthVCKMInput s 


VCKMInputs_thetal2 
VCKMInputs_theta23 
VCKMInputs_thetal3 
VCKMInputs_deltal3 


VCKM 


BlockVCKM(n) 
LengthVCKM 


VCKM_Q 

VCKM_thetal2 
VCKM_theta23 
VCKM_thetal3 
VCKM_deltal3 



Table 8: Preprocessor variables defined in SLHA.h to access the slhadata array (cont'd). 
2.2 PDG particle identifiers 

PDG . h defines the human-readable versions of the PDG codes listed in Table ^2 These are 
needed e.g. to access the decay information. At run time, the subroutine SLHAPDGName can 
be used to translate a PDG code into a particle name (see Sect. 13.12)) . 
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Block name 


Array and length 


Members 




BlockMSS2In (nl 


M^CTn M^Cfni no nl 


31,52 = 1... 3, g = 1...5 




LengthMSS2In 




7 — 
v — 




msl2in 


BlockMSL2In(n) 


MSL2In_MSL2(5i ,^ 2 ) 




MSS2In MSS2foi no 11 




LengthMSL2In 


MSL2In_MSL2Flat(i) 





MSS2In_MSS2Flat (z, 1) 


mse2in 


BlockMSE2In(n) 


MSE2In_MSE2 (31, 3 2 ) 





MSS2In_MSS2 (31 ,32 » 2) 




LengthMSE2In 


MSE2In_MSE2Flat (i) 


z 

-E- 


MSS2In MSS2Flat(? 21 

1 IkJ kJ ^1 _L XX 1 IkJ kJ ^1 X J.CLU \ v j / 


MSQ2IN 


BlockMSQ2In(n) 


MSQ2In_MSQ2(pi,5 2 ) 




MSS2In_MSS2 (31 ,32 , 3) 




LengthMSQ2In 


MSQ2In_MSQ2Flat (z) 




MSS2In_MSS2Flat (z , 3) 


MSU2IN 


BlockMSU2In(n) 


MSU2In_MSU2(pi,5 2 ) 




MSS9Tn MSSPfni no 41 

1 IkJ kJ -Lxx 1 IkJ kJ /L Vyi f Jj 2 > ' 




LengthMSU2In 


MSU2In_MSU2Flat (z) 





MSS2In_MSS2Flat (z , 4) 


MSD2IN 


BlockMSD2In(n) 


MSD2In_MSD2(^i,5 2 ) 




MSS2In MSS2(ni no 51 

1 IkJ kJ ^1 _L XX 1 IkJkJZj VyJ^ i \j £ ) ' 




LengthMSD2In 


MSD2In_MSD2Flat(i) 




MSS2In MSS2Flat(? 51 

llUU^iX XX 1 1 kJ kJ X _i_ CX O \ f y / 




BlockMSS2(n1 


l v Jooz_lJ 




= 1 ... 




LengthMSS2 


M^Qfni no nl 


9i 


32 = 1 - - - , ^ — 1 ... 5 






nooz_noozr idL \L , (jj 


i - 


= 1 ... 9, g = 1 . . . 


msl2 


BlockMSL2(n) 


mqt n 


— 


MSS2_Q CH 




LenrthMSL2 


MCI O MCJOfn, no^l 


— 










— 


TV/TO O O TV/TO O OT'H _ j_ / „" j \ 

MSS2_MSS2Flat {t , 1) 


mse2 


BlockMSE2(n1 

X-^ _1_ V ^ ill ikj'Ljij \ 1 V / 


MQT70 n 








T pntrthMSF? 


rioi_jZ_rioi-jZ vyi s y 2 / 




rlooZ_rlooZ ^31 ) 32 > 






MSE2 MSE2Flat (?'1 




MSS^FI at f? 21 

lliJiJ^. 1 lkJkJZ.1 _I_CL U ^6 j ^1 / 


MSQ2 


BlockMSQ2(n) 


MSQ2_Q 




1 ikJ kJ ^1 V'-'/ 




LengthMSQ2 


MSQ2 MSQ2((/i ,Qi) 




1 luuz 1 luuz vyl j y2 > 






MSQ2_MSQ2Flat (z) 




MSS2_MSS2Flat (i , 3) 


MSU2 


BlockMSU2(ra) 


MSU2_Q 




1 iooz 1^ v 1/ 




LengthMSU2 


MSU2_MSU2(5i,p 2 ) 




MSS2 MSS2C<7i no 41 






MSU2_MSU2Flat (z) 





MSS2_MSS2Flat (i , 4) 


msd2 


BlockMSD2(ra) 


MSD2_Q 


= 


MSS2_Q(5) 




LengthMSD2 


MSD2_MSD2(5i,^ 2 ) 





MSS2_MSS2 (31 ,32 , 5) 






MSD2_MSD2Flat(i) 


= 


MSS2_MSS2Flat(i,5) 




BlockASfMix(n) 


ASf Mix UASffsi so t) 


Sl 


s 2 = 1...6, t= 1...4 




LengthASfMix 


ASf Mix JJASf Flat (z,« 


i = 


= 1...36, t = 1...4 


SNMIX 


BlockSnMix(ra) 


SnMix JJASf (si,s 2 ) 




ASfMixJJASf (si,s 2) l) 




LengthSnMix 


SnMix_UASf Flat (z) 




ASfMix_UASfFlat(i,l) 


SLMIX 


BlockSlMix(n) 


SlMix_UASf (si,s 2 ) 




ASfMixJJASf (si,s 2) 2) 




LengthSlMix 


SIMix JJASf Flat (i) 




ASf Mix_UASf Flat (i , 2) 


USQMIX 


BlockUSqMix(n) 


USqMix JJASf (si,s 2 ) 




ASfMixJJASf (si,s 2 ,3) 




LengthUSqMix 


USqMix JJASf Flat (z) 




ASf Mix JJASf Flat (i , 3) 


DSQMIX 


BlockDSqMix(n) 


DSqMix_UASf (s 1 ,s 2 ) 




ASfMixJJASf (si , s 2 , 4) 




LengthDSqMix 


DSqMix JJASf Flat (z) 




ASf Mix JJASf Flat (i, 4) 



Table 9: Preprocessor variables defined in SLHA.h to access the slhadata array (cont'd). 
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Block name 


Array and length 


Members 


CVHMIX 


BlockCVHMix(n) 
LengthCVHMix 


CVHMix_UH(^i,/i 2 ) h!,h 2 = 1...4 
CVHMix_UHFlat (0 i = 1 ... 16 


NMNMIX 


BlockNMNMix(n) 
LengthNMNMix 


NMNMix_ZNeu(ni,n 2 ) ni,n 2 = 1...5 
NMNMix_ZNeuFlat (i) i = 1 ... 25 


NMHMIX 


BlockNMHMix(n) 
LengthNMHMix 


NMUMix_UH(/ii,/i 2 ) /ti,/i 2 = 1...3 
NMUMix_UHFlat (0 i = 1 . . . 9 


NMAMIX 


BlockNMAMix 
LengthNMAMix 


NMAMix_UA(/ii,/i 2 ) /ti, /t 2 = 1 . . . 3 
NMAMix_UAFlat (0 i = 1 . . . 9 



Table 10: Preprocessor variables defined in SLHA.h to access the slhadata array (cont'd). 



fermions 


sfermions 


PDG_nu_e 
PDG_electron 
PDG_up 
PDG_down 


PDG_snu_el PDG_snu_e2 
PDG_selectronl PDG_selectron2 
PDG_supl PDG_sup2 
PDG_sdownl PDG_sdown2 


PDG_nu_mu 
PDG_muon 
PDG_charm 
PDG_strange 


PDG_snu_mul PDG_snu_mu2 
PDG_smuonl PDG_smuon2 
PDG_scharml PDG_scharm2 
PDG_sstrangel PDG_sstrange2 


PDG_nu_tau 
PDG_tau 
PDG_top 
PDG_bottom 


PDG_snu_taul PDG_snu_tau2 
PDG_staul PDG_stau2 
PDG_stopl PDG_stop2 
PDG_sbottoml PDG_sbottom2 



bosons 




;inos 


PDG_hO 


PDG. 


.neutralinol 


PDG_HH 


PDG. 


_neutralino2 


PDG_A0 


PDG. 


_neutralino3 


PDG_Hp 


PDG. 


_neutralino4 


PDG_H3 


PDG. 


_neutralino5 


PDG_A2 


PDG. 


.charginol 


PDG_photon 


PDG. 


_chargino2 


PDG_Z 


PDG. 


.gluino 


PDG_W 


PDG. 


.gravitino 


PDG_gluon 






PDG_graviton 







Table 11: The PDG codes defined in PDG.h. 
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3 Routines provided by the SLHA library 



The file SLHA.h must be included in every subroutine or function that uses SLHALib 
routines. It contains the necessary preprocessor definitions as well as external declarations 
for the SLHALib routines. 

The basic data structure is the double complex array slhadata of length nslhadata. 
These names are hard-coded into the preprocessor definitions and may not be changed by 
the user. As a corollary, only one instance of the slhadata structure can be used in any 
one routine. This poses no serious limitation for most applications, however. 

3.1 SLHAClear 

subroutine SLHAClear (slhadata) 
double complex slhadata(nslhadata) 

This subroutine sets all data in the slhadata array given as argument to the value invalid 
(defined in SLHA.h). It is important that this is done before using slhadata, or else any 
kind of junk that happens to be in the memory occupied by slhadata will later on be 
interpreted as valid data. 

3.2 SLHARead 

subroutine SLHARead (error , slhadata, filename, abort) 

integer error, abort 

double complex slhadata(nslhadata) 

character* (*) filename 

This subroutine reads the data in SLHA format from filename into the slhadata array. 
If the specified file cannot be opened, the function issues an error message and returns 
error = 1. The abort flag governs what happens when superfluous text is read, i.e. text 
that cannot be interpreted as SLHA data. If abort is 0, a warning is printed and reading 
continues. Otherwise, reading stops at the offending line and error = 2 is returned. 
SLHARead implicitly calls SLHAClear to clear the slhadata array before reading the file. 

The blocks spinfo and dcinfo are largely ignored when reading the file, as they are for 
human information only. Only the maximum of all message codes is kept in the Severity 
member of the block. Since the message codes increase with severity, this indicates the 
overall reliability of the corresponding data (spectrum or decay information). For example, 
if the Severity member is 4 (real errors), the Accord advises not to use the corresponding 
data. See also Sect. 13.41 
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3.3 SLHAWrite 



subroutine SLHAWrite (error , slhadata, filename) 
integer error 

double complex slhadata(nslhadata) 
character* (*) filename 

This subroutine writes the data in slhadata to filename. 

3.4 SLHAInfo 

subroutine SLHAInfo (slhablock, code, text) 
double complex slhablock (*) 
integer code 
character* (*) text 

This subroutine adds a message to one of the informational blocks, SPINFO or DCINFO. 
The block is most conveniently addressed through the Block. . . macros, for example 

call SLHAInfo (BlockSPInfo (1) , 4, "Error in computation") 

Allowed codes are 

• 1 = program name, 

• 2 = program version, 

• 3 = warning message, 

• 4 = error message. 

Messages are truncated at 80 characters. 

3.5 SLHANewDecay 

integer function SLHANewDecay (slhadata, width, parent_id) 
double complex slhadata(nslhadata) 
double precision width 
integer parent_id 

This function initiates the setting of decay information for the particle specified by the 
parent_id PDG code, whose total decay width is given by width. The return value 
is an integer index which is needed to subsequently add individual decay modes with 
SLHAAddDecay. If the fixed-length array slhadata becomes full, a warning is printed and 
zero is returned. If a decay of the given particle is already present in slhadata, it is first 
removed. 
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3.6 SLHAFindDecay 



integer function SLHAFindDecay(slhadata, parent_id) 
double complex slhadata(nslhadata) 
integer parent_id 

This function also initiates the setting of decay information. Unlike SLHANewDecay, it 
requires that the decay of the parent_id particle exist and reshuffles the decay information 
in slhadata such that new channels can be added to this decay. If no decay matching 
parent_id is found, the return value is 0, otherwise it is the index needed to add decay 
modes with SLHAAddDecay. 

3.7 SLHAAddDecay 

subroutine SLHAAddDecay (slhadata, br, decay, 
& nchildren, childl_id, child2_id, child3_id, child4_id) 
double complex slhadata(nslhadata) 
double precision br 
integer decay 

integer nchildren, childl_id, child2_id, child3_id, child4_id 

This subroutine adds the decay mode 

(parent.id) — > childl.id child2_id child3_id child4_id 

to the decay section previously initiated by SLHANewDecay or SLHAFindDecay. decay is the 
index obtained from the latter (which also set the parent_id) and child?7,_id are the PDG 
codes of the final-state particles. The branching ratio is given in br. If the fixed-length 
array slhadata becomes full, a warning is printed and decay is set to zero. 

If decay is zero, an overflow of slhadata in an earlier invocation is silently assumed 
and no action is performed. It is therefore sufficient to check for overflow only once, after 
setting all decay modes (unless, of course, one needs to pinpoint the exact location of the 
overflow) . 

As with SLHAGetDecay (see Sect. I3.8|) . only the first nchildren of the childn_id are 
actually accessed and Fortran allows to omit the remaining ones in the invocation. 

3.8 SLHAGetDecay 

double precision function SLHAGetDecay(slhadata, parent_id, 
& nchildren, childl_id, child2_id, child3_id, child4_id) 
double complex slhadataO) 
integer parent_id 

integer nchildren, childl_id, child2_id, child3_id, child4_id 
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This function extracts the decay 

parent_id — > childl_id child2_id child3_id child4_id 

from the slhadata array, or the value invalid (defined in SLHA.h) if no such decay can 
be found. The parent and child particles are given by their PDG identifiers (see Sect. 12. 2|) . 
The return value is the total decay width if nchildren = 0, otherwise the branching ratio 
of the specified channel. 

Note that only the first nchildren of the childri_id are actually accessed and Fortran 
allows to omit the remaining ones in the invocation (a strict syntax checker might issue a 
warning, though). Thus, for instance, 

Zbb = SLHAGetDecay (slhadata, PDG_Z, 2, PDG_bottom, -PDG.bottom) 

is a perfectly legitimate way to extract the Z — > bb decay. 

3.9 SLHADecayTable 

integer function SLHADecayTable (slhadata, parent_id, 
& width, id, maxparticles , maxchannels) 
double complex slhadata(nslhadata) 
integer parent_id, maxparticles, maxchannels 
double precision width (maxchannels) 
integer id(0 :maxparticles , maxchannels) 

This function stores all decay channels for the particle identified by parent_id in the 
arrays id and width. Unlike SLHAGetDecay, one does not need to know the exact decay 
mode in order to extract information. The value for parent_id serves as a wildcard 
and transfers the entire decay table contained in slhadata. SLHADecayTable returns the 
number of channels found. The two arrays can be read out rather straightforwardly: 
For each channel c, 

• n = id(0,c) gives the number of participating particles, i.e. the number of decay 
products plus one. 

• The PDG code of the decaying particle is in id(l,c). 

• The PDG codes of the decay products are in id(2,c). . . id(n,c). 

• If n = 1, width (c) contains the decaying particle's total width in GeV. 

• If n > 1, width (c) contains the branching ratio for the given decay. 
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3.10 SLHAExist 



integer function SLHAExist (slhablock, length) 
double complex slhablock (*) 
integer length 

This function tests whether a given SLHA block is not entirely empty. It returns 2 if the 
block has at least one complex member, 1 if the block has at least one real member (i.e. all 
imaginary parts zero), and if the block has no valid members at all. The SLHA blocks are 
most conveniently accessed using the Block. . . and Length. . . definitions (see Sect.HJ), 
e.g. 

if( SLHAExist (BlockMass (1) , LengthMass) .ne. ) ... 

3.11 SLHAValid 

logical function SLHAValid(slhablock, length) 
double complex slhablock (*) 
integer length 

This function tests whether a given SLHA block consists entirely of valid data, i.e. it 
returns . FALSE . if at least one member of the block is invalid. The SLHA blocks are most 
conveniently accessed using the Block. . . and Length. . . definitions (see Sect. EJ, e.g. 

if( SLHAValid(BlockNMix(l) , LengthNMix) ) ... 

3.12 SLHAPDGName 

subroutine SLHAPDGName (code, name) 
integer code 
character* (PDGLen) name 

This subroutine translates a PDG code into a particle name. The sign of the PDG code is 
ignored, hence the same name is returned for a particle and its antiparticle. The maximum 
length of the name, PDGLen, is defined in PDG.h. 

3.13 Incompatible Changes 

Two incompatible changes in the interface were necessary with respect to the SLHALib 1 
[3], largely due to the fact that the SLHA2 allows complex entries: 

• slhadata is now a double complex, not a double precision array. 
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• The SLHAExist function has become an integer function, as it now distinguishes three 
possible scenarios: no valid entries, only real entries, and complex entries. 

• The SLHAWrite subroutine no longer has arguments for program name and version. 
Such informational messages can now be added with the much more general subrou- 
tine SLHAInf o. 



4 Examples 

Consider the following example program, which just copies one SLHA file to another: 

program copy_slha_f ile 
implicit none 

#include "SLHA.h" 

integer error 

double complex slhadata(nslhadata) 

call SLHARead (error, slhadata, "inf ile . slha" , 0) 
if( error .ne. ) stop "Read error" 

call SLHAInfo(BlockSPInfo(l) , 1, "My Test Program") 
call SLHAInf o(BlockSPInfo(l) , 2, "1.0") 

call SLHAWrite (error , slhadata, "outf ile . slha") 

if( error .ne. ) stop "Write error" 

end 

Already in this simple program a couple of things can be seen: 

• the file SLHA.h must be included in every function or subroutine that uses the SLHA 
routines and this must be done using the preprocessor #include (not Fortran's 
include), thus the program file should have the extension .F (capital F). 

• slhadata must be declared as a double complex array of length nslhadata. 

• One should not continue with processing if a non-zero error flag is returned. 

A more sensible application would add something to the slhadata before writing them 
out again. The next little program pretends to compute the fermionic Z decays (by calling 
a hypothetical subroutine MyCalculation) and adds them to slhadata: 
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program compute_decays 
implicit none 



#include "SLHA.h" 
#include "PDG.h" 

integer error, decay, t, g 

double complex slhadata(nslhadata) 

double precision total_width, br(4,3) 

integer ferm_id(4,3) 

data ferm_id / 
& PDG_nu_e, PDG.electron, PDG.up, PDG.down, 
& PDG_nu_mu, PDG.muon, PDG.charm, PDG.strange, 
& PDG_nu_tau, PDG_tau, PDG.top, PDG.bottom / 

call SLHARead (error, slhadata, "inf ile . slha" , 0) 
if( error .ne. ) stop "Read error" 

* compute the decays with parameters taken from the slhadata: 
call MyCalculation(SMInputs_MZ, MinPar.TB, 
& total_width, br) 

decay = SLHANewDecay (slhadata, total.width, PDG_Z) 
do g = 1, 3 
do t = 1, 4 

call SLHAAddDecay (slhadata, br(t,g), decay, 
& 2, f erm_id(t ,g) , -f erm_id(t ,g) ) 

enddo 
enddo 

call SLHAInfo(BlockDCInfo(l) , 1, "My Decay Calculator") 
call SLHAInfo(BlockDCInfo(l) , 2, "3.1415") 

call SLHAWrite (error , slhadata, "outf ile . slha") 

if( error .ne. ) stop "Write error" 

end 

Demonstrated here is the access of SLHA data (SMInputs_MZ, MinPar_TB) and the setting 
of decay information. 
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5 Building and Compiling 



The SLHA library package can be downloaded as a gzipped tar archive from the Web site 
http://www.feynarts.de/slha. After unpacking the archive, change into the directory 
SLHALib-2 . and type 

. /configure 
make 

Some simple demonstration programs (in the demo subdirectory) are built together with 
the library libSLHA.a. 

Compiling a program that uses the SLHA library is in principle equally straightforward. 
The only tricky thing is that one has to relax Fortran's 72-column limit. This is because 
even lines perfectly within the 72-column range may become longer after the preprocessor's 
substitutions. While essentially every Fortran compiler offers such an option, the name 
is quite different. A glance at the man page should suffice to find out. Here are a few 
common choices: 



Compiler 


Platform/OS 


Option name 


g77 


any 


-f f ixed-line-length-none 


pgf77 


Linux x86 


-Mextend 


ifort 


Linux x86 


-extend_source 


f77 


Tru64 Alpha 


-extend_source 


f77 


SunOS, Solaris 


-e 


fort77 


HP-UX 


+es 



To compile and link your program, add this option and -Ipath -Lpath -1SLHA to the 
compiler command line, where path is the location of the SLHA library, e.g. 

pgf77 -Mextend -I$H0ME/SLHALib-2 . myprogram.F -L$H0ME/SLHALib-2 . -1SLHA 

All externally visible symbols of the SLHA library start with the prefix SLHA and should 
thus pretty much avoid symbol conflicts. 

6 Summary 

The SLHA library presented here provides simple functions to read and write files in SLHA 
format. Data are kept in a single double complex array and accessed through preprocessor 
variables. The library is written in native Fortran 77 and is easy to build. The source code 
is openly available at http://www.feynarts.de/slha and is distributed under the GNU 
Library General Public License. 

The author welcomes any kind of feedback, in particular bug and performance reports, 
at hahn@feynarts.de. 
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